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Abstract 

An adaptive controls method for instability suppression in gas turbine engine combustors has been developed 
and successfully tested with a realistic aircraft engine combustor rig. This testing was part of a program that 
demonstrated, for the first time, successful active combustor instability control in an aircraft gas turbine engine-like 
environment. The controls method is called Adaptive Sliding Phasor Averaged Control. Testing of the control 
method has been conducted in an experimental rig with different configurations designed to simulate combustors 
with instabilities of about 530 and 315 Hz. Results demonstrate the effectiveness of this method in suppressing 
combustor instabilities. In addition, a dramatic improvement in suppression of the instability was achieved by 
focusing control on the second harmonic of the instability. This is believed to be due to a phenomena discovered and 
reported earlier, the so called Intra-Harmonic Coupling. These results may have implications for future research in 
combustor instability control. 


I. Introduction 

Lean-burning, low emission combustors are being investigated for aircraft gas turbine engines. Lean combustion 
is shown to be advantageous for reducing Nitrogen Oxides (NOx) emissions and enhancing turbine temperature 
distribution and efficiency; but is also prone to thermo-acoustic instabilities. These instabilities are typically the 
result of the coupling of the fluctuating heat release 1 (of the combustion process) with the lightly damped acoustics 
of the combustion chamber. The exact mechanisms involved in this coupling are not well understood and different 
hypotheses exist as to its precise nature. 2 

Combustor instability suppression presents a challenging problem for controls design due primarily to the large 
dead-time phase delay (of many hundreds of degrees or more) and noise in the combustion process. Besides large 
phase delay and noise, there are other characteristics of combustor instabilities such as amplitude modulations and 
net random phase walks, which could play an important role in the control design.”' 3 

A number of research efforts have attempted suppression of the thermo-acoustic instability through active 
control. 4 The goal of these active control efforts was to reduce the energy concentrated at the instability frequency 
and to reduce the overall amplitude of the combustor pressure oscillations. Some active control concepts involved 
speaker actuation 5 " 11 and others involved fuel modulation. 12 " 16 Fuel modulation is more applicable to aircraft gas 
turbine engines. Several techniques have shown some success in suppressing the frequency spectra of low frequency 
instabilities, with less success in suppressing the time domain pressure fluctuations. Almost all of the research in the 
area of active combustor instability control has made use of low-pressure, gaseous fueled combustors. One 
exception is the work by Barooah, et al. 17 that attempted to reduce pressure oscillations with the low frequency 
configuration of the same liquid fueled combustor rig described later in this article. Even though the technique 
applied had some success in reducing the instability amplitude at the instability frequency, the instability energy was 
transferred to the frequency side lobes (“peak splitting”), and the overall amplitude of the combustor pressure 
oscillations was minimally reduced. 

The Adaptive Sliding Phasor Averaged Control (ASPAC) methodology, developed by NASA and described in 
this article, has successfully demonstrated instability suppression for multiple combustor rig configurations. The 
ASPAC method was first applied to control a High Frequency Rig Configuration (HFRC) during tests conducted in 
2002 at United Technologies Research Center (UTRC). These early results, presented in reference 18, demonstrated 
the effectiveness of the methodology in reducing the combustor pressure oscillations at the instability frequency. 
This testing, which included concurrent testing of an alternate NASA control method, represented the first known 
demonstration of actively controlled, high frequency, thermo-acoustic instabilities in a realistic aircraft gas turbine 
engine combustor rig (as reported in ref. 19 and described in refs. 18, 20, and 21. In June 2004, control was 
attempted on a Low Frequency Rig Configuration (LFRC) that was assembled at NASA GRC to further test and 
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verify the control methodology. As will be shown, the LFRC exhibited a very coherent instability compared to the 
HFRC. (The actual engine instability, on which the F1FRC design was based, exhibited instability coherence 
somewhere between the LFRC and the F1FRC.) The ASPAC controller was also successful in suppressing this 
LFRC’s coherent instability. Finally, due to a previously reported coupling mechanism 18,22 involving the harmonics 
of the instability, control was then focused on the second harmonic of the instability. This approach produced 
superior results. 

This article describes the development of the ASPAC controller and the results of testing the controller with a 
combustor instability rig. The article is organized as follows. A description of the combustor rig is given that 
includes comparisons of the instabilities of the F1FRC, the LFRC, and that of a developmental aircraft gas turbine 
engine. This is followed by a description of the combustor instability behavioral model and the development of the 
control methodology. Then test results are presented for the F1FRC and the LFRC, including harmonic control and 
its implications, followed by a hypothesis that attempts to explain the phenomena observed during testing. The 
article concludes with areas of future research interests. 


II. Experimental Hardware Description 

A. Combustor Instability Rig 

In order to focus control development toward realistic combustion instabilities in aeronautics, a combustor rig 
that replicates an aircraft gas turbine engine combustor instability was designed and fabricated. 23,24 This single 
nozzle research combustor rig was developed in partnership with Pratt & Whitney and UTRC. A schematic of the 
test rig configuration is shown in figure 1 . The sample problem selected for this rig is an axial combustor instability 
that was observed during the development of a high-performance aircraft gas-turbine engine. The frequency of the 
observed instability in the developmental engine was about 525 Flz and the magnitude of the pressure oscillations 
was sufficient to cause unacceptable vibratory stresses in the turbine. 

It was previously demonstrated that the combustor rig successfully replicates the axial instability observed in the 
engine. The single-nozzle combustor rig operates at engine pressure and temperature conditions and has many of the 
complexities of the actual engine combustor. These include the same air-blast, 2-stage fuel nozzle and swirler used 
in the engine; dilution/cooling; and an effusion-cooled liner. The fuel used in the combustor rig is JetA. A venturi 
pre-diffuser and dump diffuser simulate the engine compressor exit (combustor inlet) conditions. A choked nozzle at 
the combustor exit simulates the engine turbine stator. The nearly choked venturi and choked exit nozzle, in addition 
to simulating operation inside the engine, provide acoustic isolation from the rest of the facility. 
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Figure 2. — Comparison of engine (mid-power) and 
baseline combustor rig (HFRC) pressure amplitude. 
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Figure 3. — Pressure amplitude spectra and time 
domain behavior in extended rig configuration 
(LFRC). 


Conditions corresponding to a mid-power engine condition were chosen for evaluation (T3=700 °F, P3=200 
psia, fuel-air ratio = 0.03). Water-cooled, nitrogen purged dynamic pressure transducers (PCB 12x series) were used 
to measure the pressure oscillations inside the combustor rig. Test results (fig. 2) established the existence of a 
combustor instability at about 533 Hz — approximately the same frequency as the developmental engine. The test rig 
apparatus for the 533 Hz HFRC is shown in figure 1. The LFRC is obtained by removing the pre-diffuser section 
shown in figure 1 and placing it before the two 14 wave spool sections. This has the effect of elongating the 
combustor by approximately 19 in. and produces an instability at 315 Hz. 

For comparison, the pressure spectrum in the actual engine and in the single-nozzle combustor rig at comparable 
operating conditions for the HFRC and LFRC is shown in figures 2 and 3 respectively. This comparison shows that, 
in terms of instability frequency and amplitude, the HFRC closely simulates the actual engine instability. In terms of 
instability coherence for the fundamental (i.e. instability magnitude with respect to noise), however, the engine more 
closely resembles the behavior of the LFRC. These two rig configurations provide essentially two control problems. 
The HFRC provides a high frequency, low signal-to-noise instability problem, and the LFRC provides a lower- 
frequency, high amplitude, high signal-to-noise problem. 

The first experimental testing with the HFRC took place at UTRC, and the more recent control testing with the 
LFRC took place at NASA GRC. 

B. Fuel Modulation Valve 

For the open-loop fuel modulation testing and closed-loop instability control testing, a high-frequency fuel 
modulation valve developed by Georgia Institute of Technology was used. The valve is based on a magnetostrictive 
material actuator. The valve was typically located approximately 18 in. (45.7cm) from the fuel injector tip. The 
open-loop combustor pressure response to fuel valve modulation for two different runs is shown in figure 4. The 
response was somewhat nonlinear in that larger valve commands produced a disproportionately larger response in 
combustor pressure. However, the response of the fuel modulation system was determined empirically to be 
satisfactory for open-loop testing and for closed-loop instability control. The dynamics of the aforementioned fuel 
modulation sub-systems were not investigated further. This is, however, deemed to be an important area for future 
research. The fuel valve control system also had a mean flow control loop which kept the mean flow constant while 
the fuel flow was perturbed. This ensured that the results seen were due to fuel modulation and not to a change in 
fuel flow and/or fuel-air ratio. The correct operation of the mean flow control was verified during test and from post- 
test data. 
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Figure 4. — Normalized combustor pressure response to perturbation 
of the high frequency fuel modulation valve. 


III. Combustor Instability Simulation 

The ASPAC control diagram is shown in figure 5. A reduced-order behavioral model is used to represent the 
combustor instability. This model consists of the flame dynamics, acoustics, and a saturation non-linearity. The 
instability is self-excited, that is, it requires no input via fuel injection modulation. The respective transfer functions 
for the flame (Gp), acoustics (Ga), and non-linearity (NL) are: 


G F = 


Kpcop 
s + co F 


(1) 


G a = 


K A (Q 2 a 

s 2 + 2 g A a> A + 


(2) 


NL = tw\i(p) (3) 

where p is the acoustic pressure output from the transfer function in Eq. (2). 

In this simulation, the flame dynamics in Eq. (1) are represented with a first order transfer function with 
a>P »l.lo>^and a gain Kp =1.2. The acoustics transfer function dynamics in Eq. (2) correspond to an instability 
frequency of ~ 533 Hz for the HFRC or ~ 315 Hz for the LFRC. With g A = 0.0707 and K A = -0.25 , the coupled 
response between the flame and acoustics becomes oscillatory and unstable (growing with time). The purpose of the 
non-linearity (NL) is to limit the amplitude of the instability. These choices of the transfer function parameters, with 
the addition of any noise magnitude other than zero, allow for self sustained oscillations that closely approximate the 
combustor instability from time response measurements. 4 ' 12 Figure 6 shows a simulation of a self-excited instability 
with some injected white noise added to represent the process noise in the combustor. The response is also 
numerically unstable in the sense that the response depends on the damping g A and the sampling time T s 
( T s = 0.0001 sec), however, the results are deterministic and repeatable. 

The dynamics of the fuel valve, feed lines, injector, and heat release or flame dynamics due to the fuel 
modulation (Gy in Fig. 5), have been simulated with a second order transfer function as 
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( 4 ) 


Gy = 


r v 


S 2 + 2 gyCOyS + CDy 


with K v =1.0 and its natural frequency ©^comfortably greater than ©4 so that it will not attenuate the controlled 
fuel modulations. g v is chosen to provide a damped response. 
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Figure 5. — Combustor instability control block diagram. 
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Figure 6. — (a) Self excited instability time response, 
(b) Waveform of fully developed instability. 


NASA/TM- 


-2008-215202 


5 






For the filter shown in figure 5, two identical band pass filters are placed in series. The purpose of these filters is: 
(a) to pass the measured instability for control; and (b) when running the control with a simulated combustor, to 
simulate the large dead-time phase shift found in the combustion process. This phase shift amounts to approximately 
760° for the LFRC and is significantly larger for the F1FRC. The state space model of the band pass filter is: 


x = Ax + Bu 
y = Cx 


( 5 ) 
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For this simulation, a\ HFRC = 0.8 , a2 HFRC = 3.58, b HFRC = 800 and a\ LFRC = 0.4 , a 2 LFRC = 1.69, b LFRC = 400. The 
Bode plot of the band pass filter for the LFRC is shown in figure 7. 


IV. ASPAC Algorithm and Modifications 

The ASPAC algorithm is described in detail in references 18 and 25. In principle, the algorithm calculates a 
restricted control phase region in a stationary frame of reference, which is favorable for instability suppression, as 
shown in figure 8. First, the combustor instability pressure is sensed using a band-pass filter in order to isolate the 
instability frequency. The pass band of the filter is based on some a priori knowledge of the approximate instability 
frequency and its width covers the expected instability frequency drift within the operating range of an actual engine. 
The sensed pressure is phase shifted, with a phase that slides back and forth inside a restricted control region. This 
phase shifted pressure is used to generate a command to the fuel valve which produces a fuel modulation opposing the 



Frequency, Hz 


Figure 7. — Bode plot of band pass filter. 
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instability pressure. The dashed line region (so called effective stability region) within the shaded area in figure 8 is 
the stable region in which the phase of the fuel modulation relative to the instability pressure is such that the power 
of the combustor pressure oscillations is reduced. The combustor pressure resulting from heat addition due to the 
fuel modulation, vectorially adds to the combustor pressure due to the instability to provide the overall dynamic 
combustor pressure. 

One of the difficulties in this control problem is that the effective stability region continues to shrink as the 
instability is being suppressed, and eventually disappears, rendering the system unstable. This causes the stable 
region to open up again and the process is repeated. Combustor noise plays a major role in this process in terms of 
how far the instability can be suppressed. 2 At higher suppression levels, the noise increasingly drives the process 
causing the instability phase to vary abruptly. 

The pressure sensing and application of the control phase shift occurs at a rate of 10 KHz, which amounts to a 
phase resolution of approximately 18° for the HFRC instability. The controller calculates a new phase shift at a rate 
of 40 Hz. This slower control rate is due to the large dead-time phase shift, which necessitates an appropriately 
reduced controller bandwidth. Also, this controller optionally employs Discontinuous Exponential Gain Modulation 
Control 1 * 25 (DEGMC). In this control mode, the gain is modulated on/off with an exponential decay in order to 
counteract the effective proportional gain variability produced by the large dead time phase delay of the plant. 
DEGMC was improvised for this particular application and it is not a standardized control gain modulation 
technique. In addition, controller parameter adaptation 18 22 is employed in order to adapt some of the key parameters 
of the controller. 

The control design was later modified to also include harmonic control. 22 The reason for this modification was 
evidence of the so called intra-harmonic coupling mechanism, discovered by analyzing test results. 218 Figures 9 and 10 
show this analysis. First, figure 9 shows that the instability pressure grows even though the modulated fuel command is 
opposing it. This instability growth persists beyond the time due to the large dead-time phase delay of the plant, and 
suggests another phenomenon at work. Figure 10 shows a strong coherence between the fundamental and the second 
harmonic of the instability in the frequency region of about 200 to 400 Hz suggesting a coherent energy transfer 
between the fundamental and the harmonic. Both these results have been produced by certain manipulations of the test 
data described briefly here and more fully in references 2 and 18. Figure 9 shows the combustor rig pressure and the 
controlled fuel modulation command for the HFRC. Both signals are filtered to show only the instability and the 
respective control command. The fuel control command is shown opposing the instability rather than the combustor 
pressure due to the fuel modulation because the latter cannot be measured separately. However, this combustor 
pressure will be proportional to its control command, displaced by a fixed phase. As a result, figure 9 was produced by 
sliding the time series of the control waveform relative to the overall combustor rig pressure so that the phases of the 
two signals will be opposing at the regions where the instability pressure is decreasing (this much is known). This 
manipulation effectively changes the control signal to an approximation of its corresponding proportional combustor 
pressure. Even though the controlled signal phase starts opposing the combustion pressure soon after the pressure 
oscillations begin to grow again, the pressure oscillations continue to grow for a while. It was suspected, based on 
results like figure 9, that the fundamental frequency and the harmonics are somehow coupled as inputs/outputs through 
an unknown mechanism. Figure 10 shows this analysis for the LFRC. One cannot compute the coherence of signals of 
different frequencies. Therefore, the original combustor pressure signal was band-pass filtered to isolate the 
fundamental instability mode. Then the original signal was also re-sampled in order to halve its frequency components. 
The re-sampled signal was again filtered to isolate the “fundamental mode” except this time, the “fundamental mode” 

Pressure 
from fuel 
modulation 


Boundary 
of effective 
stability 
region—' 


Figure 8. — ASPAC phasor diagram in 
a stationary frame of reference. 
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is actually the original second harmonic of the instability. Finally, the coherence of these two resulting signals was 
computed. This effectively gives the coherence between the fundamental and the second harmonic as seen in 
figure 10. As shown, the coherence between the second harmonic and fundamental (in the frequency band around 
300 Hz) is high, around one. In figure 10, the coherences of some of the higher order harmonics are shown as well 
because the band-pass filters do not completely get rid of the higher order harmonics. The coherence suggested that 
additional suppression of the fundamental instability frequency may be possible by also simultaneously controlling 
the harmonics. 

As a result of this analysis, the controller structure was modified by adding a 2 nd harmonic control and 3 rd 
harmonic control to take advantage of the intra-harmonic coupling effect. These modifications, which also include 
the parameter adaptation portion of the control algorithm, are shown in figure 11. The change with regard to 
harmonic control allows selecting the controller structure so as to focus on the fundamental and/or higher harmonics. 
The limitation of this control structure is the frequency that can be influenced by the fuel valve, that is, control can 
only be focused on the harmonic frequencies if the fuel valve has bandwidth/authority at those frequencies. 


V. Controller Test Results 

This section describes simulations and test results obtained using the ASPAC control methodology. Simulations 
were conducted initially to help develop the algorithm and to predict the controller behavior to the extent possible 
before the start of testing. Two separate tests were then conducted to validate the methodology. The first set of tests 
was carried out on the engine-like HFRC at UTRC, in June and September 2002. The second test was performed on 


the LFRC at NASA GRC, in June 2004. 



Figure 9. — Filtered combustor instability pressure versus 
filtered and phase shifted fuel valve control command. 



Figure 10. — Instability harmonics coherence plot. 
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Figure 1 1 . — ASPAC combustor instability control block diagram for harmonic control. 
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Figure 12. — Simulation of combustor instability 
(a) control response and (b) power spectral 
density after instability suppression. 


A. Simulations 

Figure 12(a) shows the time response and figure 12(b) shows the power spectral density of the combustor pressure 
with the combustor instability simulated using Eqs. (1) to (5). The instability is allowed to grow before control is applied at 
about 0.6 sec. With no a priori knowledge of the phase required for suppression, the controller very briefly causes the 
instability to grow. The controller then quickly identifies the proper control phase to apply and suppresses the instability to 
reduce the instability pressure oscillations. Figure 12(b) shows the power spectral density when the controlled instability 
has been suppressed to about -50 dB (near the noise floor for this simulation), where the uncontrolled instability 
magnitude is about 0.6 psi or -4.4 dB. Brief control simulations were also conducted against a sectored 1-D physics based 
model of the combustor process 26 in order to further validate the controls design. 

B. HFRC Test Results 

The objective of the F1FRC control testing was to demonstrate control of combustor instability under conditions similar 
to those experienced in an actual aircraft gas turbine engine (e.g., similar in terms of the operating conditions, instability 
frequency, and instability amplitude, see fig. 2). As described earlier, the band pass filter used in the simulation also 
simulates the large phase shift experienced in the combustor instability process. Flowever, for rig testing, the filter design 
was changed to remove most of the phase shift shown in figure 7. Recall that the reasons for selecting the filter design for 
the simulation were to help isolate the instability frequency, and also to simulate the process phase shift. As a result, a 
reduced phase shift filter with a 450 to 700 Flz pass band was selected for testing: 
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Through simulation and experimentation, it was found that best control results are achieved with a somewhat wide 
band pass band filter for both the HFRC and the LFRC (i.e., about 250 FIz and 100 Hz pass band respectively). 

Test results using the aforementioned control algorithm for the HFRC are shown in figure 13. Combustor pressure was 
sensed approximately 1 in. downstream from the fuel injector face using the previously described PCB dynamic pressure 
transducers. The fuel flow was modulated with the Georgia Tech fuel valve. The control algorithm was implemented on a 
dSpace modular control system. In both the June and September 2002 tests, the controller successfully suppressed the 
amplitude spectrum of the pressure oscillations at the instability frequency to near the noise floor. In terms of time domain 
reductions, the instability amplitude in the HFRC is small compared to the overall wideband noise (~l/7 ,h the noise 
amplitude) and no apparent reductions were visible in the time domain. As seen in the upper plot of figure 13, a lower 
frequency coupling (~5 Hz) existed in the combustor rig during the June 2002 testing. It was observed that occasional 
coupling between the combustor pressure and the facility open-loop fuel feed pressure was causing the facility fuel 
pressure to oscillate at this low frequency. This was remedied for the second test through filtering and tuning of the mean 
fuel flow controller. Figure 14 shows test results with and without DEGMC. Some noticeable improvement 
(approximately 25% reduction) was observed by utilizing the DEGMC portion of the algorithm. Figure 13(b) shows the 
results from the second HFRC test with the fuel feed system control and DEGMC improvements. 
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Figure 13. — Amplitude spectral density of uncontrolled 
versus controlled instability; (a) June and 
(b) Sept. 2002 testing. 
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Figure 14. — Instability suppression with and without 
discontinuous exponential gain modulation control. 


C. LFRC Test Results 

The objective of the LFRC control testing was to demonstrate the control algorithm on a combustor rig with 
operating conditions similar to an actual aircraft gas turbine engine, but with an instability which is more coherent 
(i.e., less damped) than that of the HFRC. This more closely approximates the instability coherence of the actual 
engine (see figs. 2 and 3). Initially, during testing with the LFRC, control was focused on just the fundamental 
frequency of the instability ( — 315 FIz). Through this test, good control of approximately 90% suppression of the 
amplitude spectrum at the instability frequency was obtained with a corresponding 35% reduction in time domain 
peak pressure oscillations and a 70% reduction in RMS pressure (see fig. 15 compared to fig. 3). In addition to the 
combustor instability pressure (PLA1C1), figure 15 also shows the controller fuel valve command, in volts 
(REFAcmd). As shown earlier, the fuel modulation command causes a modulation pressure in the combustor of 
approximately 1 psi/volt up to about 350 Liz. For this test, approximately 20% fuel pressure modulation was 
required for instability suppression. 

During the same test run, controller action was focused on the second harmonic. The purpose of this test was to 
see if the instability at the fundamental frequency responds to control of the second harmonic (as discussed earlier), 
in order to attempt additional instability suppression. It soon became apparent that the instability fundamental not 
only responded to the second harmonic control, but the suppression of the fundamental was more drastic than first 
expected. As shown in figure 16 the peak instability pressure amplitude is 0.7 psi. Comparing this result to the 
uncontrolled instability of figure 3 there is ~95% suppression in amplitude at the instability frequency, a ~75% 
reduction in time domain peak pressure oscillations, and an 80% reduction in RMS pressure. The corresponding fuel 
valve command for this controlled combustor pressure is also shown in figure 16. From earlier frequency sweep 
testing, it had been observed that a 2 volt peak-to-peak fuel valve excitation at 600 FIz was barely registering a 
noticeable pressure response (i.e., pressure above the noise floor). Therefore, the relatively little control authority at 
the second harmonic seems to be all that was needed to suppress the instability. The band pass filter used to focus 
control effort on the second harmonic is wide enough to allow a small amount of effort at the fundamental. 
Flowever, the higher frequency control action is what seems to provide the improved results over the fundamental- 
only control. 


NASA/TM— 2008-215202 


11 



Run 423, Date 040527, Point 019 


Run 423, Date 040527, Point 055 



Figure 15.— Amplitude spectra and time domain of Figure 16.— Amplitude spectra and time domain of 

controlled instability pressure (PLA1C1_psi) and controlled instability pressure (PLA1C1_psi) and 

control command (REFAcmd). control command (REFAcmd) with harmonic 

control applied — Case 1. 


The results from another control test also focusing on the second harmonic are shown in figure 17. For this test, 
parameter adaptations were enabled in the controller. These adaptations included primarily the controller gains and 
the parameters for the exponential gain modulation control. This result is even better in terms of the peak amplitude 
spectra than the result shown in figure 16 (0.39 versus 0.7 psi peak, 97% suppression as compared with fig. 4). In 
both tests the second and higher harmonics of the instability have been almost completely eliminated. Some peak 
splitting is evident in this latter test. Depending on the particular controller design, evidence of peak splitting is more 
pronounced at increased levels of suppression where combustor noise becomes more dominant and net random 
phase walks become more abrupt. 2 It is thought that this phenomenon could translate to frequency changes 
impacting energies at the instability side lobes. A comparison of the uncontrolled (as in fig. 3) versus the controlled 
instability for this last test is shown in figure 18. 

Testing focusing on the third harmonic of the instability was also carried out during the course of this test run. 
However, based on this testing and from frequency sweep tests, very little fuel modulation authority of the fuel 
valve exists at the third harmonic of ~945 Hz (no visible modulation authority is observable from frequency sweeps 
tests at these high frequencies). As a result the third harmonic control testing was only able to reduce the instability 
by a few psi, down to about 8 psi peak pressure (results are not shown) at the instability frequency (as compared to 
>13 psi for the uncontrolled case). 

The significance of second harmonic control is that better instability suppression seems to be achievable through 
this control, as the results presented here seem to support. The drawback is that harmonic control may require a 
higher fuel actuator bandwidth, and possibly a penalty on the fuel actuator life due to higher cycling rates. Both 
problems could possibly be mitigated by the seemingly little actuator authority needed to suppress the instability at 
the second harmonic, with correspondingly small magnitude dithering of the valve and small actuation rates. It could 
be possible that similar results may also be obtained by exercising sub-harmonic control, or control at a significantly 
lower frequency than the instability. The next section will offer some discussions along these lines. Interaction with 
other instability modes (e.g., tangential modes) would be an interesting research subject as well, but was not 
investigated as part of this effort. 
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Run 423, Date 040527, Point 056 



Frequency, Hz Time, sec 

Figure 17. — Amplitude spectra and time domain of 
controlled instability pressure (PLA1C1_psi) 
and control command (REFAcmd) with 
harmonic control applied — Case 2. 



Frequency, Hz 


Figure 18. — Amplitude spectra of uncontrolled versus 
controlled instability for the LFRC. 


VI. Discussion of Results 

While the primary result of the effort reported in this article is the successful suppression of combustion 
instability in an aircraft gas turbine engine type combustor, a discussion of how this was accomplished (focusing 
control effort on the second harmonic) and why this was effective seems important as well. In this section a 
hypothesis is offered to explain the mechanism that causes the phenomena of “Intra-Harmonic Coupling” in 
combustor instabilities. Also, discussed is the possible role of impedances in the root cause of combustor 
instabilities. These discussions are based on the previously discussed instability control test results and also on other 
phenomena observed during testing. 

A. Cause of Intra-Harmonic Coupling — Hypothesis 

In addition to closed-loop controls testing, additional test points were conducted where the combustor pressure 
was perturbed at various discrete frequencies, and also where semi-continuous frequency sweeps were carried out. 
This was done mainly to characterize the fuel valve modulation authority. However, in the process of conducting 
these tests, some interesting phenomena were observed. Applying a discrete modulation of 2 or 3 volts peak, at a 
lower frequency (in the range of 200 Hz), initially, caused a barely visible response in the pressure amplitude 
spectra. But after about 10 minutes of running time, this modulation would rapidly grow to greater than 2 psi in the 
combustor pressure as if some reorganization of the combustor oscillatory waveform was taking place (fig. 19). The 
growth of this discrete modulation seemed to interact with the instability, reducing its peak somewhat. Further, 
discrete excitations at frequencies within a few Hz on either side of the instability frequency of approximately 
315 Hz had an even more pronounced effect. While the peak at the excitation frequency would grow significantly, 
the peak at the instability frequency would be reduced dramatically (fig. 20). This latter effect can be explained by 
what can be called instability entrainment, whereby the system is susceptible to oscillations and will exhibit 
oscillatory behavior at frequencies near its resonance when excited at these frequencies; thereby increasing the 
system’s energy at the driven frequency and at the same time reducing the energy at the resonant frequency. 

Applying a discrete fuel modulation at frequencies below 200 Hz showed another interesting phenomena. For a 
100 Hz modulation frequency (but not at the 155 Hz half-frequency of the instability), a strong interference with the 
instability was observed (fig. 21). 

For all these open- loop modulation cases, it was observed that discrete frequency excitations seemed to create 
their own sub-harmonics, as well as harmonic frequencies. Overall, this combustor seemed to be very conducive to 
dynamic coupling between harmonics and also between discrete frequencies. 
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Run 424, Date 040601, Point 019 


Run 424, Date 040601 , Point 024 



0 500 1000 0.0 0.1 0.2 0 500 1000 0.0 0.1 0.2 


Frequency, Hz Time, sec Frequency, Hz Time, sec 

Figure 19. — Amplitude spectra and time history of Figure 20. — Amplitude spectra and time history of 

combustor pressure (PLA1C1_psi) after 10 min combustor pressure (PLA1C1_psi) with 289 Hz 

with 200 Hz Fuel perturbation command fuel perturbation command applied (REFAread). 

applied (REFAread). 


Run 425, Date 040602, Point 023 



Frequency, Hz Time, sec 

Figure 21 . — Amplitude spectra and time history of 
combustor pressure (PLA1C1_psi) with 100 Hz 
fuel perturbation command applied (REFAread). 

All these behaviors, involving couplings between harmonics and between different frequencies, could be caused 
by the manifestation of different dynamic processes. But it might make more sense to consider that, for the most 
part, the same underlying mechanism is involved and can explain all these phenomena. In general these coupling 
behaviors could be caused by Wave Coupling Asymmetry’ (WCA) as was described in some sense in reference 27 
that deals with Oceanography. In terms of the combustion process, this involves the coupling of pressure waves of 
different frequencies. By WCA, the meaning here is that when two or more waves interfere, the interference is non- 
symmetrical, causing a net suppression or amplification of the individual waves. 

Let’s consider for a moment energies of different frequencies, like two pressure waves, existing inside a 
confined domain like a combustor chamber and oriented in the same spatial direction. Even though these waves are 
of different frequencies, they can still interfere with one another. Such interference will manifest as cycle-by-cycle 
suppression or amplification of each wave, depending whether their phases are such that they promote amplification 
or suppression at a given instant in time. However, over the course of a combined periodic cycle (defined by the 
least common multiple of the two frequencies), the net change in suppression and amplification of each wave will be 
zero. Therefore, in an amplitude spectra density plot, the two waves will seem unaffected by each other. This can be 
expected to be the case, unless there is asymmetry (uneven coupling) in either the direction of suppression or 
amplification. 
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In the case of the combustor pressure waves; when the phases of the two waves are such that amplification is 
supported, this action is likely opposed by damping in the form of friction (like skin friction with the walls of the 
chamber). This friction damping could be nonlinear, with damping increasing as the amplitude of the wave 
increases. On the other hand, when the phases of the two waves are such that it promotes suppression in some 
portion of the overall periodic cycle, the suppression will be uninhibited by damping. In such a case there will be 
WCA or unbalance in suppression and amplification which would favor increased suppression with time. However, 
as suppression increases, this asymmetry may no longer favor suppression. This is because as the wave amplitude 
diminishes and becomes comparable to the amplitude of the combustor noise, the noise increasingly drives the 
process thus preventing further suppression. At this point WCA then favors amplification. At some point a balance 
is reached, where the two sided WCA reaches an equilibrium point. In an amplitude spectral density plot, the two 
waves reach a balance between suppression and amplification, where likely the higher amplitude wave looses some 
net energy and the lower amplitude wave gains some net energy. This could explain what was observed in figures 19 
to 21. 

The coupling of the discrete frequency excitation, which results in the growth of the discrete frequency (after 
some ten minutes of running time on the test rig) and the simultaneous suppression effect of the instability, cannot, 
however, be explained by the WCA alone. Perhaps a Minimum Energy ; System (MES) effect is also involved with 
these coupling phenomena. In such a case the relative phase of the two waveforms could also be evolving with time 
to cause suppression of the higher amplitude waveform and some amplification of the lower amplitude waveform, in 
order for the system to settle into a MES state. 

If this hypothesis is correct, the question remains why controlling the second harmonic has such an immediate 
and major impact in suppressing the fundamental mode. The answer could be related to the periodicity of the 
combined wave cycle. In the case of the fundamental and its harmonics, this combined cycle has the same 
periodicity as the fundamental. This means that an action to suppress the fundamental can be applied as frequent as 
once every cycle of the fundamental. However, this would mean that applying a control action to suppress the 
fundamental via the second harmonic, would inevitably result in applying an equal action to amplify the 
fundamental, within the same cycle. But as was discussed before, coupling asymmetry due to damping, favors 
suppression. And because of the high frequency of control action when working with the harmonics, large levels of 
instability suppression are possible down near the noise floor where balance of the two sided WCA is reached. This 
integer multiple of periodicity between the fundamental and its harmonics in this combustor environment that favors 
wave interference based on the WCA, could be the cause of the intra-harmonic coupling mechanism discussed 
before. 

The observed interaction between harmonics begs the question: Can sub-harmonic control could be used 
effectively to suppress the fundamental mode of the instability? In other words, the confined domain of a combustor 
chamber may promote coupling of energy between pressure or thermo-acoustic waves of different frequencies. But 
the coupling between waves with integer multiple periodicities may be more efficient, and therefore stronger. 
Interfering with this coupling as was done by focusing control on the second harmonic, seems to result in improved 
instability suppression, in fact, altogether causing the elimination of higher order harmonics which seem to behave 
as interdependent coherent structures. The possibility of sub-harmonic control (focusing control on frequencies that 
have integer harmonics at the instability fundamental) is hinted at in the result shown in figure 21. Unfortunately, 
testing of this hypothesis could not be investigated due to facility schedule limitations. 

With a more pronounced or coherent instability where the damping is less, the level at which equilibrium in the 
WCA is reached, will likely be higher. Meaning that with more coherent instabilities, the WCA will still favor 
suppression but to a lesser degree. Based on this, the combustor noise alone (given the large dead time phase delay 
in the process), would not seem to be the only dominant factor in combustor instability suppression of coherent 
instabilities. The other limiting factor is likely the instability damping or the instability coherence. 

B. Possible Root Causes of Combustor Instabilities 

Some of the unexplained experimental observations concerning the open- and closed-loop response to fuel 
modulation bring up the more fundamental topic of the root cause of combustor instabilities. The predominant 
explanation in the literature is that combustor instabilities are caused by thermo-acoustic coupling, that is, coupling 
of the heat release with the acoustics of the combustor chamber. Based on this understanding, some experimenters 
have attempted to adjust the relative phase between the heat release and the acoustics by adjusting the fuel transport 
delay time, in order to passively suppress the instability. This approach has had mixed results, 28 ' 29 which suggests 
that the relative phase is not the only mechanism responsible for the presence of an instability. In reference 2 it was 
discussed that even though the thermo-acoustics are responsible for providing the resonance that promotes pressure 
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oscillations in the combustor chamber, the exact nature of the instability, like magnitude and frequency, would seem 
to be the result of the fluidic impedances involved in this process (i.e., source/load type impedance interactions). 
Impedances play a major role in the dynamic behavior of energy transfer systems, especially in systems that are not 
decoupled through choking of orifices. Therefore, a complete understanding or characterization of the instability 
behavior can not ignore the impedances involved in this process. This is becoming an area of investigation as 
reported in the literature, although mainly for premixed, gaseous fueled gas turbines. 30,31 Gaining a more in depth 
knowledge of the role that impedances play in this process, especially for liquid fueled, aircraft gas turbine engine 
type combustors, is of interest to the authors as a topic of future research as it may allow for more effective control 
designs. 

The discussions in this section suggest several hypotheses concerning the nature of the mechanisms responsible 
for the dynamic couplings of combustor instabilities. The aim of this discussion is to promote further investigation 
through experiments and analysis that will either confirm or negate such hypotheses. In any case more knowledge 
will be gained about the behavior and control of combustor instabilities. 

VII. Concluding Remarks 

This article reports on test results of an adaptive control method for suppression of combustion instabilities. 
Results are presented from two different test runs. The first set of runs utilized a high frequency configuration of a 
combustor rig which emulates an aircraft gas turbine engine combustor instability. The second set of test runs 
utilized a low frequency configuration of the same rig. The test results validate the effectiveness of this adaptive 
control method in suppressing combustor instabilities in liquid-fueled gas turbines. A certain characteristic called 
“Intra-Harmonic Coupling,” discovered in previous analysis was exploited to produce even higher suppression 
levels in the control of the instability. This was accomplished by focusing control at the second harmonic of the 
instability. A hypothesis has been offered to explain the mechanisms involved in intra-harmonic coupling, which is 
attributed to a two sided wave coupling asymmetry influenced by the combined wave periodicity, and also possibly 
a minimum energy systems effect. Areas of interest for future research are: to apply the demonstrated phase shifting 
control method to advanced, low-emissions combustors; to develop actuation systems with more predictable 
dynamic response; and to develop a better understanding of the harmonic coupling effects and other instability 
mechanisms in order to support the development of even more effective control strategies (e.g., impedance -based 
control). Ultimately the goal of these combustor instability control efforts is to enable NOx emissions reductions 
through the application of advanced control approaches to ultra-low-emissions, multi -point, lean direct injection 
combustors. 
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